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ABSTRACT 

Most magnetic bearing control schemes use a bias current with 
a superimposed control current to linearize the relationship between 
the control current and the force it delivers. With the existence of the 
bias current, even in no load conditions, there is always some power 
consumption. In aerospace applications, power consumption becomes 
an important concern. In response to this concern, an alternative mag- 
netic bearing control method, called Adaptive Variable Bias Control 
(AVBC), has been developed and its performance examined. The 
AVBC operates primarily as a proportional -derivative controller with 
a relatively slow, bias current dependent, time- varying gain. The 
AVBC is shown to reduce electrical power loss, be nominally stable, 
and provide control performance similar to conventional bias con- 
trol. Analytical, computer simulation, and experimental results are 
presented in this paper. 

INTRODUCTION 

Magnetic bearings are electromechanical devices that use mag- 
netic forces to completely levitate a rotor without physical contact. 
The forces they produce depend nonlinearly upon the current flow- 
ing in the magnet coils and the distance between the rotor and the 
stator poles. In conventional operation of a magnetic bearing, a 
method known as the bias method is used to linearize the relation- 
ship between the current supplied to and the force delivered by the 
bearing. These devices cannot operate without closed loop electronic 
control because they are open loop unstable due to its attractive- 
only nature. Usually classical control techniques such as proportional- 
derivative (PD) or proportional-integral-derivative (PID) feedback 
control are used. 

For many space applications, conservation of power is a high 
priority issue. The research presented here focuses on the develop- 
ment, application, and effectiveness of Adaptive Variable Bias Con- 
trol (AVBC) to reduce the power required to operate a magnetic bear- 
ing, consequently resulting in power-saving, while preserving stability 
and actuator frequency response. Reduced bias operation of mag- 
netic bearings would not only provide electrical power-saving but 
would provide many other synergistic benefits, such as reduced power 
requirements, reduced heat transfer problems, and reduced rotor 
heating due to eddy currents and hysterises. This method is not pur- 
ported to be a new method for rotordynamic control but rather a 
modification to the operating mode of a magnetic bearing to make it 
a more efficient control actuator. Analytical, simulation, and experi- 
mental results are presented. 

Low or no bias research has been done by the primary author as 
part of his Ph.D. dissertation (see Johnson, 1995), involving zero 
bias and pseudo-zero bias control. Zero bias work has also been done 
by Brown and Grodsinsky, 1991, Ishida, 1990, and Johnson, et al., 


1991. Some early work was done by Cunningham, et al., 1978, and 
Sperry Right Systems, 1979, that involved limited real-time bias 
adjustment. Extensive analysis of power loss (hysterises, eddy cur- 
rent, windage, configuration, etc.) in magnetic bearing systems has 
been done by Kasarda, et al., 1997. Recent work reported involves 
an advanced control method applied to low bias systems (Knospe 
and Yang, 1997). 


ACTIVE MAGNETIC BEARING OPERATION 
Magnetic Bearing Force Equations 

A typical magnetic bearing comprises a set of radially posi- 
tioned electromagnets positioned in opposing pairs around a lami- 
nated magnetic bearing journal (See Fig. 1). For a magnetic bearing 
with four electromagnets there is one opposing pair for each perpen- 
dicular axis. Each electromagnet consists of a laminated core and 
one or more coil windings. In the experimental configuration used 
in this research the poles are uniformly positioned 45 degrees apart 
(See Fig. 2). 

The force produced by a single two pole electromagnet can be 
shown to be given by the following equation where / is the total cur- 
rent in the magnet coils, Xg is the gap distance, is the permeabil- 
ity of free space, A is the pole face area, and N is the number of coil 
turns: 


^p/UV 2 / 2 


( 1 ) 


The force in Eq. (1) is attractive and increases as the gap 
decreases. This attractive force produces an unstable system for an 
open loop magnetic bearing configuration. The net force, F n , pro- 
duced by an opposing pair of identical two-pole electromagnets on a 
single axis is the sum of the forces produced by each electromagnet; 
taking account of the sign convention (see Fig. 1 ), the net force equa- 
tion is given as follows: 
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lj is the current in magnet j, and Xgj is the gap distance for magnet j. 

Previously, the equations were given for the static force pro- 
duced by a magnetic bearing along a single axis. That force is af- 
fected dynamically by the rate limit at which current changes in the 
coils, called the current slew rate limit, which is dependent on the 
voltage limit, V max , of the power supply and the coil inductance, L. 
An approximate expression for the slew rate limit is equal to 
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Conventional Bias Current Operation 

Conventional magnetic bearing practice linearizes the force- 
current relationship by using a bias current along with a superim- 
posed control current. In order to obtain a linear relation between 
force and control current we express the currents 7j and I2 in the 
opposing coils as: 


l l =I b + , c 

l2=Ib~Ic 


(3) 


where is the bias current and l c is the perturbation or control cur- 
rent. Next we express each gap in terms of the nominal gap Xq 
and the deviation X from this nominal value caused by rotor 
displacement: 


In the defining equations below, capitalized variables are dimen- 
sional and the lower case variables generally are nondimensional. 
Additionally, subscripts 1 and 2 refer to the right and left magnets, 
respectively. 
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I r€ f is the reference, or nondimensionalizing, current, / max is the 
maximum allowable coil current, and V max is the maximum power 
supply voltage. 

We define a bias critical frequency, Q cr as the frequency at 
which voltage saturation first occurs for the conventional bias method 
with the bias current equal to the reference current: 


X gl =X o-X 


Xg2 =X$+X 
Combining Eqs. (2) to (4) yields: 


(4) 
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The frequency ratio is defined as follows, where Q is the com- 
mand frequency: 


Fx 
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A linear relationship between force and control current and position 
can be obtained by assuming X « Xq and l c « /^, which can be 
shown to result in the following equation: 


4 Z/f ATI. 

F n = — 5^-Xh = K x X + KiI c (6) 

*0 X 0 

K x is the bias (negative) stiffness and X, is the current stiffness. 

This linearized force equation is generally used for linear con- 
trol development even though the assumptions that were made are 
frequently violated in practice. For proportional-derivative feedback 
control the following relation is used. 


«0 = ~- (10) 

il cr 

The nondimensional time, T, whose units are cycles, is given as 
follows: 


t = Qt 


(11) 


The instantaneous power, P f -, is defined as follows: 


I) = l 2 R 

The total energy, E { , is defined as follows: 
^rad^ 

E t = ji 2 Rda rad t 
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( 12 ) 


(13) 


I c =-(K p X + K d x) (7) 

where Kp is the proportional control gain and K ^ is the derivative 
control gain. 

ADAPTIVE VARIABLE BIAS CONTROL 
Preliminary Control Design 

To simplify the equations describing the controls, the impor- 
tant variables will be transformed into nondimensional variables. 
The nondimensional factors will be used for both AVBC and con- 
ventional bias control methods but are derived from parameters and 
variables that are predominantly associated with the bias method. 
They are also based on the force and power electronic equations 
presented earlier. 


T is the period the total energy is being evaluated over and the sub- 
script rad denotes radians. 

The total nondimensional, e t , energy is defined as follows: 


e t = 


The average nondimensional power (total energy per cycle), p t will 
be defined by the next equation. This variable will be used to evalu- 
ate the power used by each method. 
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Conventional Bias (CB) 

The conventional bias method produces a linear operating ac- 
tuator which enables use of many linear control techniques. The 
means by which this is accomplished was outlined earlier. The bias 
current is set to one-half the designed maximum allowable current 
to ensure that the control current never exceeds the bias current. If 
the command current magnitude exceeds the bias current, linearity 
remains intact but, due to the attractive-only behavior of each mag- 
net, the resultant force would be counteractive to the desired actua- 
tion action. The bias current is responsible for the major portion of 
power dissipated when operating the actuator using this method. The 
nondimensional bias force, f c ^ is: 

Fx^ 

fcb=-T- = *c 06) 

zr , 


Adaptive Variable Bias Control (AVBC) 

A new proposed operating method for magnetic bearings is the 
Adaptive Variable Bias Control method (AVBC), based on the Con- 
ventional Bias (CB) method. The primary idea behind this control 
method is that by varying the bias current adaptively as needed, the 
power requirement can be reduced while producing the desired force 
and maintaining the linearization afforded by using bias. The main 
static influence in changing the bias current is to change the bearing 
static stiffness. Dynamically, there are two effects: the control laws 
depend on bias, and as the bias changes, the force slew rate limit 
changes. Since there are so many possibilities as to how to effec- 
tively vary the bias, the use of adaptive control is introduced. Adap- 
tive control entails monitoring the operation of the magnetic bearing 
and determining from that information the best setting for the bias 
current for those conditions. Adaptive variable bias appears to be 
efficient and flexible enough to consider for a variety of magnetic 
bearing applications. The application of adaptive control here is a 
loose application of the self-tuning type of nonlinear adaptive con- 
troller. Here the bias current is a system parameter that can be changed 
by using information about how the system is operating. The “esti- 
mation’* of the bias is not based on typical parameter estimation but 
it is determined by an online estimate based on the commanded cur- 
rent. We require that the bias must be equal to the amplitude of the 
command current. Thus, a simple adaptation law is obtained in that 
we find an average command current amplitude over some time span 
and set the bias equal to this average. 

The operating control law can best be described as follows: 

'l=|'c| + , c- '2=|'c|-'c ( 17 > 

and the resulting nondimensional adaptive variable bias force, f av fr 

is: 

favb =4 M‘c (18) 

EXPERIMENTAL SET-UP 

In order to demonstrate the AVBC method on an actual physi- 
cal system, experiments were performed. The set-up of the 
experimental magnetic bearing system is comprised of several inter- 
connected electrical and mechanical components: a magnetic bear- 
ing, the power amplifiers, an IBM PC AT compatible computer, the 
A/D and D/A boards, and a Rotor Kit. The Rotor Kit, which is 
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manufactured by the Bently Nevada Corporation, is a precision model 
of a high-speed rotating machine that is designed for laboratory dem- 
onstration of radial shaft vibration. It can operate at speeds up to 
12 000 rpm. Its general use is for demonstrating experimental rotor 
dynamic phenomena but it is being used here with magnetic bear- 
ings to demonstrate the control methods presented in this work. The 
magnetic bearing is installed on the rotor kit at the end of the rotor 
away from the motor. The magnetic bearing is an 8-pole magnetic 
bearing (see Fig. 3) and was fabricated at the NASA Lewis Research 
Center. The design specifications are given in Fig. 2. The force out- 
put of the magnetic bearing was characterized by the single force 
parameter, Z, determined experimentally to be approximately 
5.76e-4 lbf-(in./A)2 w hj c h w in be the Z value used when calculating 
force from experimental current. Important parameters are listed in 
Table 1. 

RESULTS 
Conventional Bias 

The conventional bias method is governed by the equations 
developed earlier. The simulations were run at four different force 
command levels with the command frequency equal to the bias criti- 
cal frequency (frequency at which saturation occurs for peak force 
equaling maximum force output) of 80 Hz. Figure 4 is a plot of 
several cycles of the resulting nondimensional currents due to vari- 
ous desired nondimensional forces for each electromagnet of the 
horizontal axis. For the current with the largest amplitude (dotted 
line) there is a noticeable linear portion as the positive peak ampli- 
tude is approached. This linear portion is due to voltage saturation. 
In Fig. 5, a plot of the force produced and the command current 
based on the commanded force is shown. A comparison of the com- 
mand current in Fig. 5 and the resultant currents in Fig. 4 show that 
the desired output is achieved unless voltage saturation occurs; then 
the output is less than desired. The comparison of the output force 
and command current show that they are only different by a constant 
factor. Voltage saturation occurs only for the force with amplitude 4. 
These simulation results are virtually the same as those expected 
from the analytical description. 

The experimental results for CB operation, with 2 A bias cur- 
rent and a particular command amplitude at the bias critical frequency 
of 80 Hz, is shown in Fig. 6. Figure 7 shows the current, calculated 
force, and measured force for the bearing operating at low 
amplitude and frequency force command. This verifies the force pro- 
ducing ability under CB operation of the bearing with respect to 
analytical prediction. 

Adaptive Variable Bias 

The Adaptive Variable Bias Control (AVBC) method is gov- 
erned by the equations developed earlier. The simulation of the AVBC 
method was run at the same four force commands as with the con- 
ventional bias. In Fig. 8 is a plot of several cycles of the resulting 
nondimensional currents for these four cases. Note that the form of 
the currents changes differently with commanded force than for the 
CB method (see Fig. 4). Also notice that the dc component of the 
current is approximately equal to the amplitude of the ac compo- 
nent. Hence for every force amplitude the minimum value of the 
currents is always near zero. Additionally, as with conventional bias 
for/= 4, there is an appearance of linear form due to voltage satura- 
tion. In Fig. 9 a plot of the force produced and the command current 
based on the commanded force is shown. A comparison of the com- 
mand current in Fig. 9 and the resultant currents in Fig. 8 show that, 
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in general, the desired output is achieved; but if voltage saturation 
occurs the output is less than desired. The comparison of the output 
force and command current show that they are only different by a 
constant factor, except when voltage saturation significantly affects 
the force at an amplitude of 4. The simulation results yielded the 
variable bias operational characteristics as expected from the ana- 
lytical description. 

The experimental results for AVBC operation are shown in 
Figs. 10 to 12. Figure 10 shows the resulting current for a particular 
command amplitude at the bias critical frequency of 80 Hz. Fig- 
ure 1 1 shows the current, calculated force, and measured force for 
the bearing operating at low amplitude and frequency force com- 
mand. This verifies that the force produced under AVBC operation 
matches the analytical prediction. Figure 12 shows how this method 
changes the bias current according to the control current, indicating 
very stable adaptation features. 

Force and Power Consumption Evaluation 

In Fig. 1 3, there are plots of analytical, experimental and simulat- 
ed power versus force data for the AVBC and the CB method. For 
the CB method, the data reveals a parabolic relationship between 
force and power with a nonzero power value at zero force, attributed 
to the ever-present bias current. The analytical solution matches well 
with simulation data. The simulation data point that is farthest to the 
right is off from the analytical results because the simulation in- 
cludes the resistive effects in the current dynamics which adds to the 
voltage causing voltage saturation earlier than for inductive voltage 
only and for high current values it becomes more influential. The 
experimental results match the other results well. For the AVBC 
method, there is a linear relationship between force and power with 
a zero power value at zero force. The analytical solution matches 
well with simulation data, except for the point mentioned before. 
The experimental results again match the other results well. 

DISCUSSION 

The results of the analysis of the Adaptive Variable Bias Con- 
trol (AVBC) data has been presented. The analytical results were 
based on linearized equations, the simulation results were based on 
nonlinear equations, and the experimental results were from real- 
time, physical results. These results indicate how a magnetic bear- 
ing operates using the AVBC method. The Conventional Bias (CB) 
and AVBC methods produced the commanded force for the majority 
of the range of force commanded. For high force command there 
was some voltage saturation for simulated and experimental results 
for these bias-based methods. The effectiveness of the AVBC method 
to operate with minimal power was very good, in general, compared 
to the CB method. The CB method does not do well for low force 
levels due to the invariable bias current. As the commanded force 
level increases the difference in power dissipation between each 
method begins to decrease. For /> 3.5, the CB and AVBC methods 
have comparable power levels. The results presented in this work 
shows that the AVBC has performed just as well as the CB method 
in producing force and is significantly better in using less power. 

CONCLUSION 

In conclusion, this research effort provided an analysis of AVBC, 
with a focus on operating with minimal power, and showed that it 
worked. This method was compared to the conventional operating 
method and it was determined to be viable. Further development can 


be performed by considering other adaptation schemes. The AVBC 

method proved to be a method of choice for power saving operation. 
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TABLE 1 . — POWER ELECTRONICS AND OTHER 
EXPERIMENTAL PARAMETERS 


Amplifier bandwidth 

12 kHz 

Amplifier gain 

0.5 A/V 

Power supply current limit 

10A 

Power supply voltage limit 

±35 V 

User set maximum current 

4Aa max ) 

User set maximum voltage 

lOVCV^) 

Number of coil turns, N 

236 

Cross-sectional pole area, A 

0.35 in. 2 

Nominal air gap length, x 0 

0.022 in. 

Position probe sensitivity 

200 V/in. 

Inductance of magnetic coil, L 

10 mH 

Resistance of magnetic coil, R 

0.9 Q 

Air permeability constant, jIq 

2.82e-7 lbf/A 2 

Force transducer constant 

2 Ib/mV 
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Figure 1.— Typical magnetic bearing configuration 
showing electromagnets only on one axis. 


r- 1 .650 ±. 005R 



Figure 2. — Experimental configuration of an eight pole 
magnetic bearing. 



Figure 3. — Eight-Pole Radial Magnetic Bearing imple- 
mented in Bentley- Nevada Rotor Kit in Magnetic 
Bearing Research Lab at NASA Lewis Research 
Center in Cleveland, Ohio. 


Nondimensional force command amplitude 
f=1 f = 3 




Figure 4. — Nondimensional output current for bias 
method simulation, (a) Right magnet, (b) Left magnet. 


Nondimensional force command amplitude 


f = 1 f = 3 




Figure 5. — Bias method simulation, (a) Nondimensional 
output force, (b) Command current. 
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Results 




Figure 6. — Conventional bias experimental nondimen- 
sional current for command current i c = .375. (a) Right 
magnet, (b) Left magnet. 


Nondimensional force command amplitude 


f = 1 f = 3 




Figure 8. — Nondimensional output current for adaptive 
variable bias method simulation, (a) Right magnet. 

(b) Left magnet. 




Force 


Predicted 

Measured 




Figure 7. — Conventional bias experimental dimensional 
results for command current l c = 0.0227 amps and 
ft = 8.9 Hz. (a) Right magnet, (b) Left magnet 
(c) Command current, (d) Measured and predicted 
force. 


Nondimensional force command amplitude 


f = 1 


f = 3 



Nondimensional time 



Figure 9. — Adaptive variable bias method simulation, 
(a) Nondimensional output force, (b) Command 
current. 
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Results 


1 



Nondimensional time 



Figure 10. — Adaptive variable bias experimental nondi- 
mensional current for command current i c = .375. 

(a) Right magnet, (b) Left magnet. 




Force 

Predicted 


Measured 




Figure 1 1 . — Adaptive variable bias experimental dimen- 
sional results for command current l c = 0.203 amps 
and ft = 8.9 Hz. (a) Right magnet, (b) Left magnet. 

(c) Command current, (d) Measured and predicted 
force. 



700 


Nondimensional time 



Nondimensional time 



Nondimensional time 


Figure 1 2.— Adaptive variable bias experimental currents, 
(a) Nondimensional command current, (b) Right magnet 
control current, (c) Left magnet control current showing 
adaptive characteristics. 
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results 
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bias method 
results 


— Analytical 

x Experimental 

o Simulated 


Analytical 

O Experimental 
□ Simulated 


o 

>» 

0 

1 

i 

o 

a 

« 

c 

o 

to 

c 

<D 

E 

b 

c 

o 



Nondimensional force 


Figure 13. — Nondimensional power per cycle versus non- 
dimensional force for the conventional bias method and 
adaptive variable bias method with t*> = 1 . 
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